Pulsed Raman output coupler for an atom laser 
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We theoretically study a pulsed stimulated two-photon 
Raman outcoupler for an atom laser using a full three- 
dimensional description. A finite-temperature trapped Bose- 
condensed atomic gas is treated self-consistently by the 
Hartree-Fock-Bogoliubov equations. The model is closely re- 
lated to a recent experiment on optical outcoupling [E.W. 
Hagley et al., Science 283, 1706 (1999)]. We analyze the mo- 
mentum distribution of the output atoms and show how the 
output beam may be used as a probe of the quantum state 
for the trapped atomic gas and how it could be engineered 
and controlled in a nonlinear way. 
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The generation of Bose-Einstein condensates (BECs) 
of weakly interacting atomic gases has the potential to 
open a whole new area in atom optics; coherent atom 
optics. To do this, a well-controlled source of coher- 
ent matter waves needs to be developed. Rudimen- 
tary coherent matter wave sources ("atom lasers") have 
been engineered experimentally [1,2], but there is a need 
for sustained development. In this Letter we present a 
quantum- mechanical, three-dimensional (3D), ab initio 
treatment of a pulsed Raman outcoupler for a BEC. We 
demonstrate how the output of the atom laser can be 
engineered and how it may be used as a probe of the 
quantum and thermal properties of the BEC. We use a 
well-tested self-consistent field-theoretical treatment to 
describe the finite-temperature, multi-mode BEC (the 
analogy of the cavity in a conventional laser), and well 
controlled and detailed approximations for the output 
coupling to a quasicontinuum of states. 

A fundamental difference between optical and atom 
lasers is that atoms interact while photons do not. This 
significantly complicates the theoretical analysis of atom 
lasers, as compared to optical lasers. Studies using clas- 
sical mean fields (e.g. Gross-Pitaevskii equation) disre- 
gard thermal fluctuations, decoherence, and information 
about quantum statistics. Due to the inevitable com- 
putational difficulty of a complete 3D analysis, there 
has not existed before a rigorous 3D, multi-mode, finite- 
temperature description for an outcoupler of an atom 
laser, despite a notable research activity on atom laser 
output theory in recent years [6-12]. 

We study a pulsed outcoupling which is amenable to 
a simpler theoretical description than a general continu- 
ous coupling, since it allows us to reach the Markovian 
limit and to ignore any nontrivial coherences between the 



trapped and free modes. A non- Markovian output may 
emerge in the continuous case [7,12], which, in a multi- 
mode description for the trapped gas, can lead to very 
complicated nonequilibrium dynamics. 

We start from the Hartree-Fock-Bogoliubov (HFB) 
[3,4] field theory for the Bose-condensed atoms, harmoni- 
cally trapped in internal level |c), in thermodynamic equi- 
librium. The Hamiltonian density for the trapped atomic 
gas alone reads 
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where Hem 



/(2m) + muP'T^ 12 is the center-of- 
mass Hamiltonian in an isotropic trap with frequency uj. 
The interaction strength = Anach^ /m is given in terms 
of the 5-wave scattering length ac, and /i^^^ denotes the 
chemical potential. In the HFB theory we decompose the 
condensate and the noncondensate fractions into mutu- 
ally orthogonal subspaces and expand the field operator 



(2) 



where the quasiparticle creation and annihilation oper- 
ators a}^ and an obey the standard bosonic commuta- 
tion relations [3]. The BEC wave function ^(r) satisfies 
the generalized Gross-Pitaevskii equation with the ther- 
mal population acting as an external potential [4]. We 
self-consistently solve for the mode functions Un and Vn 
and the BEC density profile (as explained in detail in 
Ref. [4]) in such a way that the Hamiltonian Hi is di- 
agonal in the Popov approximation (where two-particle 
correlations due to the anomalous average of the fluc- 
tuating field operator are neglected). The HFB-Popov 
treatment has been exceedingly successful in describing 
the properties of BECs in both isotropic traps [4] and, at 
relatively low temperatures, also in anisotropic traps [5]. 

We study the outcoupling of the trapped quantum de- 
generate Bose-Einstein gas from internal sublevel |c) = 
|^,m) to an untrapped level \o) = \g^m"). We con- 
centrate on a situation where the atoms in level \o) 
do not experience a magnetic trapping potential. As 
in the NIST experiments [2] we consider a two-photon 
output coupling via two noncopropagating lasers: Level 
|c) is optically coupled to an electronically excited state 
|e) = |e,m') by the laser field Ei with the frequency 
r^i, and state |e) is coupled to level \o) by the driving 
field E2 with the frequency (^2- The atomic transition 
frequencies for |c) \e) and |e) \o) are uji and UJ2 
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{cOoc = ^1 — ^2), respectively. The lasers are off-resonant, 
so that the population in level |e) remains small and 
spontaneous emission may be ignored. In the rotating 
wave approximation (rwa) it is useful to define the slowly 
varying field operator ipo = e*^^^"^^)^^^ fQj. IqyqI |o) in 
the Heisenberg picture and the positive frequency com- 
ponents E+ = e*^^^E+, where E+ = ejSj{rt)/2. In the 
general formalism we allow for the lasers to be focused at 
the center of the trap, although the main results in this 
paper are obtained for nonfocused lasers with a Gaus- 
sian time profile: £j{rt) oc /^^^ ^e*^-^'^. We set the 
positive z axis to be the direction of the momentum kick 
(q = ki — k2). For pulsed outcoupling with strong mo- 
mentum kick perpendicular to gravity [2], gravity is not 
expected to play as important role as in a rf output [1]. 
Nevertheless, in the diagnostics of the trapped BEG, it 
would be advantageous to cancel completely the effects 
of gravity. We assume this is done, e.g., by means of 
a pulsed magnetic field gradient or an optical potential 
during the output pulse and ignore gravity in the follow- 
ing analysis. 

We adiabatically eliminate the excited state and keep 
only the terms of first order in 1/Ai, the inverse of the 
one-photon detuning Ai = uji — Qi. With the outcou- 
pling the Hamiltonian density is Hi + W2, where 
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and Soc = ^oc — + is the two-photon detuning. 
Hem = — 7i^V^/(2m), and the light-induced level shifts 
Si and the Rabi coupling k, are defined by 
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Here the reduced dipole matrix element di {d2) for the 
atomic transition |c) \e) (|o) |e)) also contains the 
corresponding nonvanishing Clebsch-Gordan coefficient. 

By inserting Eq. (2) into the coupling terms propor- 
tional to K in Eq. (3) we immediately observe that, de- 
spite introducing the rwa, the creation of one output 
atom may be associated with the annihilation (V^Jo^n) or 
creation ('^JceJ^) of a quasiparticle in a trap [11]. Gonse- 
quently, even at T = 0, and in the weak coupling limit, 
the outcoupling of an atom can create elementary excita- 
tions in the quasiparticle vacuum demonstrating how the 
full multimode approach with quasiparticle excitations is 
unavoidable in obtaining the complete description of the 
outcoupling process. 

We consider a low density for the outcoupled atoms so 
that /i^^^ c:^ and we may ignore the term proportional 
to go in Eq. (3). Moreover, we expand the output matter 
field in the plane wave basis ipo{'^) = Sk(^l^)^k5 with 
(r|k) = exp (ik • r), where Ok is the annihilation 

operator for the mode k. We may then obtain the Heisen- 
berg equation of motion for slowly varying output mode 
operators Ok = exp [i(ek + (j^oc)t]o\^ with ek = Tik^/ (2m): 



We have defined the coupling matrix elements, analogous 
to the Franck-Gondon factors, by y^^\t) = (k|/^(rt)|/) 

andak'(t) = (k|[(52(r^)-^,,n(^)(r)/n]|k0e^(^^-^^')^ The 
coUisional interaction between the trapped and output 
atoms is treated semiclassically: ^pl^poi^li^c iplipon^^^ . 

For the Markov and Born approximation to be valid we 
ignore the depletion of the trapped modes, requiring that 



the light only weakly perturbs them, S-\ 



< 



and that we 



only couple out a small fraction of atoms, (3k <^1. In the 
first Born approximation we write an{t) c± Q^n(— cxo)e*^^^. 
For a continuous outcoupling the Markovian limit would 
be more restrictive, since then the atoms should also 
leave the interaction region (the trapped gas of size R) 
much faster than any timescale r required to form co- 
herences between the trapped and free modes, indicating 
rhq/m ^ R. The timescale r determined by the col- 
lision rate r ~ T^/ {Qoc'^^^^) or by the coupling rate l/n 
necessitates a very strong laser kick q. 

In Eq. (5) ^kk' describes the coupling between the dif- 
ferent output modes as a result of the laser focusing 
and the interactions with the trapped atoms. For short 
pulses, P <C %/ {qoc'^^^^)^ (or when Ti?q/m ^ Rgoc^^^^) 
the atomic collisions between the trapped and free atoms 
can be ignored during the laser coupling, and we may set 
goc — in ^kk' [13]. Moreover, if also the spatial vari- 
ation of the laser intensity is negligible over the atomic 
cloud (nonfocused case), the coupling between the dif- 
ferent k modes vanishes due to the orthogonality of the 
plane waves and we obtain ^kk' — ^2^(k — k^). In this 
case the solution of Eq. (5) reads: 



lim OkW/(2i7r) = J2 h-'^"" + S',, - 

t — >oo ' L (3 

n 

-aJ,Z>i*^(52 (ek + C + ^n) (^k + O- (6) 

Here the operators an and are evaluated at the ther- 
modynamic equilibrium t ^ —00 and 5'^^ = 5oc — /J^^^^ /h. 
The dynamics of the coupling matrix elements has been 
factored out: 



/ 



e^V/3V}^)(t) 



ly^iO). We 

also introduced a delta- function with a finite width F: 
^r(^) = 7r~-^/^r~^ exp (— x^/r^). In the limit of long 
pulses P ^ 00 the width is zero and Sr{x) becomes the 
Dirac delta function. 

In Eq. (6) we have ignored the light-induced level 
shifts in the delta functions by assuming a weak cou- 
pling 6i ^ cjn, 6k. The physical interpretation of Eq. (6) 
is straightforward. The expression represents the output 
mode amplitudes after the laser pulses, P ^ t. The con- 
servation of the momentum is expressed in the coupling 
matrix elements u between the modes. The integrand of 
ly is nonnegligible over some finite volume, with the char- 
acteristic length scale Rq, resulting in a momentum un- 
certainty Ak ^ I/Rq. The delta functions dictate the en- 
ergy conservation. Due to the finite duration of the pulse 
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the delta functions have acquired a finite width, inversely 
proportional to the pulse length; again a direct conse- 
quence of the Heisenberg uncertainty relation. The sign 
of the quasiparticle energy in the delta functions is differ- 
ent in the terms proportional to Un and corresponding 
to the annihilation and creation of an excitation associ- 
ated with the outcoupling of an atom. In the limit of a 
short pulse, /3^(ek + ^oc^^n)^/4 <C 1, the delta functions 
become constant, and we obtain Ok(oo) :^ ipn^^'^i'^^ (0)^ 
where tpc is evaluated in the initial state t ^ —oo. Then 
the output after the pulse is proportional to a one-particle 
correlation function of the trapped atoms: 

{oio^) oc J dVdV'e^^^-^)-^'^-'^')(V^t(r)^^(r')) . (7) 

For a homogeneous BEC the output mode density can 
measure a Fourier component of the standard first-order 
spatial coherence function. Equation (7) should be con- 
trasted to the Bragg spectroscopy which probes two- 
particle correlations of the trapped atoms [14]. 

In a more general case the output mode density is ob- 
tained from Eq. (6). If the collisions between the trapped 
and the output atoms cannot be ignored, or if the laser 
is focused, the result (6) is no longer valid, and we need 
to solve the coupled set of differential equations (5). In 
general, evaluating the output field boils down to cal- 
culating the dynamics and the coupling matrix elements 
^u^\ etc. In a spherically symmetric trap we may expand 

the mode functions Un{r) = T.nlrnCnlfh<l^nli^)'^lfh{^ 

[similarly for Vn{'^)] in terms of the radially symmet- 
ric wave functions 4>fii{r) and the spherical harmonics 
Yi^{0,ip). The symmetry is simplified in the case of 
nonfocused lasers or when the focused lasers are par- 
allel. Then the outcoupling is cylindrically symmetric 
around the direction of the momentum kick {z axis). 
Consequently, the Rabi frequencies i<i{r^ 0) are indepen- 
dent of the polar angle Lp which can be integrated analyt- 
ically: K^i^ = {k\f^{r,0)\un) = {kp,kz\K(r,e)\(t)^i), where 
we may write {p,z\kp,kz) = Aj^Pj^^z /r)Jfh{kpp)e'^-^ , 
with = x'^ -\-y'^. Here Jk denotes the kth order Bessel 
function of the first kind, the associated Legendre 
function, and Af^ complex coefficients. As a result, the 
momentum distribution of the atom laser is character- 
ized by the two components, kz and /Cp, representing the 
momentum along the laser kick and perpendicular to it. 

In the numerical calculations we consider an outcou- 
pling process with nonfocused lasers from an isotropic 
3D trap with the trapping frequency co = 27t x lOOHz by 
evaluating the operator expectation values for the output 
mode density from Eq. (6). The BEC profile, the ther- 
mal population, the excitation frequencies cj^, and the 
quasiparticle mode functions, Un and Vn^ are solved for 
10000 ^^Na atoms with the scattering length a = 2.8nm 
and the coupling strength ^ = /^(0,0) = O.OScj. This 
involves a self-consistent evaluation of about 1000 exci- 
tation frequencies and the respective mode functions [4]. 



The (2 X 2)D coupling matrix elements /^^/-^ are numeri- 
cally integrated for (n, T, m). With 10000 atoms the num- 
ber of outcoupled atoms in some examples remains small, 
but it becomes experimentally reasonable if we scale the 
results for experiments on several millions of atoms. 

At T = 60nK we obtain 7070 BEC atoms in the trap 
with the corresponding 2930 atoms in the noncondensate 
fraction. In Fig. 1 we display the momentum distribution 
of the outcoupled atoms corresponding to the different 
coupling processes proportional to <l>, Un^ and v*. We 
have chosen the detuning S'^^ = —20a;, the pulse length 
P = 0.5/cc;, and the momentum kick q = 0.1 x 47r/A, 
where A = 589nm represents the wavelength of light 
for typical optical transitions in ^^Na. The momentum 
distributions corresponding to the three different mech- 
anisms are clearly separated due to the different reso- 
nance energies. This results in the spatial separation of 
the three output fields, which could be used to probe the 
quantum properties of the trapped gas. 
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FIG. 1. The momentum distribution of the outcoupled 
atoms along the direction of the momentum kick (kz) and per- 
pendicular to it (kp) representing the coherent output from 
the BEC (left), the annihilation (middle), and the creation 
(right) of a quasiparticle in the trap. The corresponding pro- 
cesses are peaked at kzl — 5.6, 6.8, and 5.9 resulting in a 
spatial separation of the propagating output beams. 

Energy conservation, determined by Eq. (6), is opti- 
mized when 6g -h -h — hqk^^ /m-\-5'^^ = ±uJn , where k^^ 
denotes the momentum of the atoms along the laser kick 
relative to the momentum kick of the lasers, /c^ = q — kz- 
With q = {) the output is spherically symmetric. Energy 
conservation is optimized for different values of k'^ for 
the output processes proportional to <I>, or 'u*, with 
—ix)n and uon representing 'u* and respectively. The 
energy condition is in general matched at two different 
values of k'^. In the case of Fig. 1 we obtain for the BEC 
k'^ ^ -1.9// and k'^ ^ ll/l [I = (fi/mcj)^/^], represent- 
ing the atoms flying in the opposite directions from the 
trap. However, in general, the two energy maxima can 
represent very different occupation numbers due to the 
different strengths of the coupling matrix elements. In 
particular, momentum conservation for the coherent cou- 
pling from the BEC is optimized when k'^ = kp = 0, i.e., 
when eq = —5'^^- -"-^ -^^S- 1 o^^Y the first maximum is dis- 
played which is approximately three orders of magnitude 
larger than the second one. Nevertheless, it is fascinat- 
ing that, with large BECs, we still find parameter regimes 
where the output from a BEC mode is represented by two 
distinguishable fragments with an observable number of 
atoms flying in the direction opposite to the laser kick, 
indicating a coherent splitting of the BEC output. 
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Whether the output distribution is more determined 
by the optimized energy or momentum conservation de- 
pends on the spatial confinement of the trap and on the 
duration of the laser pulses /3. In Fig. 2 we show the 
momentum distribution for atoms coupled out of the 
BEC mode for three different values of P describing a 
transition from well-optimized momentum conservation 
to well-optimized energy conservation. 

As the final example we investigate the number of out- 
coupled atoms Nq. In Fig. 3 we show the number of 
coherently and incoherently outcoupled atoms as a func- 
tion of the pulse length for q = 0.1 x 47r/A and for 
different values for the detuning. We observe that Nq 
may also decrease with increasing f3. This seemingly 
counter-intuitive result can be explained by noting that 
the uncertainty of energy conservation decreases with in- 
creasing pulse length. The finite-width delta function 
in Eq. (6) then approaches the Dirac delta function and 
a nonresonant coupling becomes suppressed. Even very 
small changes of the laser detuning Sqc can result in large 
modifications of the output. Therefore, e.g., externally 
induced level shifts could possibly provide a strongly non- 
linear control of the atom laser output. 




FIG. 2. The momentum distribution of the coherently out- 
coupled atoms from the BEC for three different values of the 
pulse duration f3, with fixed Soc — — lOcj, N = 7073, and 
q = 0.05 X 47r/A. With /3a; = 0.6 (left) momentum conserva- 
tion is well-optimized representing a peak at kz — q — 2.2/1. 
For longer pulses energy conservation becomes better opti- 
mized resulting in arising second peak at kz — 2q. With 
Pco = 0.8 (middle) the two peaks exhibit equal hights. With 
Poo = 1.0 (right) energy conservation dominates. 

The essential features of the coherent output in Fig. 3 
may be understood by deriving an approximate expres- 
sion for the total number of atoms, coherently coupled 
out from the BEC, by assuming a Gaussian density pro- 
file for the BEC with the width R along the laser kick: 

No - NR^(3\exp [-x^f3\e, + Soc)^^] , (8) 

where = {qRY /[{qRY ^^{P^qY]- In the derivation of 
Eq. (8) we assumed that the width of the BEC density 
profile perpendicular to the laser kick satisfies Rp/l ^ 
PuoX' The approximate expression (8) describes well the 
functional dependence of the output. In Fig. 3 the best fit 
is obtained by choosing R ~ 0.5Rtf for S^^^ < —lOuo and 
R ~ O.GRtf for > —lOu^ where the Thomas- Fermi 
radius Rtf = l{l5Nac/iy^^. For short pulses 2{peqY < 
{qR^ ^ we obtain x — 1, and the output depends quadrat- 
ically on the pulse length A^o ^ Nn^P'^. For long 
pulses, 2{Peq)'^ {qR)'^^ the output is linear with 



p: No - NR^pqR/{V2eq)exp{-[qR{1^6oc/eq)m. 
There is also an interesting intermediate regime, where 
(qR)'^ - 2{peq)'^ and ^ 1/2. Provided that the cou- 
pling is sufficiently off-resonant, so that we simultane- 
ously have 2 ^ P\^q -\- Soc\^ the exponential factor in 
Eq. (8) becomes important, and the number of output 
atoms may rapidly decrease as a function of the pulse 
length for fixed coupling strength. This behavior results 
from energy conservation and the very narrow energy 
width of the BEC and is a consequence of the genuine 
quantum nature of the BEC, as compared to thermal 
atomic ensembles with a broad energy distribution. 
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FIG. 3. The number of coherently outcoupled atoms from 
the BEC mode (left) as a function of the pulse length p. 
Beginning from the lowest curve, the value of the detun- 
ing changes evenly from S^c/^ = —20 to —10 {cq :^ lOo;). 
The coherent coupling dominates the output close to T = 0. 
The right plot shows the coherently outcoupled atoms (solid 
lines) and the total number of outcoupled atoms at T = 60nK 
(dashed lines) for b'od^ — —16, —14, and —10. 
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